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High-statistics differential cross sections for the reactions yp — > pr\ and 7p — > pr/ have been 
measured using the CLAS at Jefferson Lab for center-of-mass energies from near threshold up to 
2.84 GeV. The rj results are the most precise to date and provide the largest energy and angular 
coverage. The rj measurements extend the energy range of the world's large-angle results by approx- 
imately 300 MeV. These new data, in particular the rj measurements, are likely to help constrain 
the analyses being performed to search for new baryon resonance states. 

PACS numbers: ll.80.Cr 11.80.Et 13.30.Eg 14.20.Gk 



I. INTRODUCTION 

Studying low-energy r\ and r( photoproduction 
presents an interesting opportunity to search for new 
baryon resonances. Since both of these mesons have 
isospin 1 = 0, the Nij and N-q' final states couple to 
N* states but not A* states. Previous experiments 
have produced precise cross section measurements for the 
IP ~ * P 7 ! reaction from threshold up to a center-of-mass 
energy, W, of approximately 2.5 GeV [3, 0, H, 111- For 
the 7p — > pr]' reaction, previous results are fairly precise 
from threshold up to W w 2.2 GeV ||. 

Studies performed on these previous experimental data 
have yielded evidence for nucleon resonance contribu- 
tions. For example, Anisovich et al. Q confirmed that 
rj photoproduction is dominated near threshold by con- 
tributions from the 5n(1535) and 5*11 (1650) states. Ev- 
idence was also found for contributions from other res- 
onances at higher energies, along with strong i-channel 
contributions in the forward direction. The previously 
published CLAS rj' results support contributions from 
several resonance states as well [5j . 

The rj results presented in this paper are more precise 
than any previous measurements and extend the large- 
angle energy range by approximately 600 MeV. They will 
provide stronger constraints on models that attempt to 
extract resonance contributions in this reaction. Our rj 
measurements extend the energy range of the world's 
large-angle results by approximately 300 MeV. Signifi- 
cant discrepancies exist between our rj results and those 
previously published by CB-ELSA [2J at higher energies 
(see Section [VIII|) . These new results for both rj and rj 
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will surely have an impact on the physics interpretation 
of these reactions. 



II. EXPERIMENTAL SETUP 

The data were obtained using the CEBAF Large Ac- 
ceptance Spectrometer (CLAS) housed in Hall B at the 
Thomas Jefferson National Accelerator Facility in New- 
port News, Virginia. A 4 GeV electron beam hitting a 
10~ 4 radiation length gold foil produced real photons via 
the bremsstrahlung process. The recoiling electrons were 
then analyzed using a dipole magnet and scintillator ho- 
doscopes in order to obtain, or "tag" , the energy of the 
photons [7f. Photons in the energy range from 20% to 
95% of the electron beam energy were tagged and thus 
measured with an energy resolution of 0.1% of the elec- 
tron beam energy. The data were analyzed in center-of- 
mass energy bins that varied in width from 10 MeV up 
to 40 MeV, depending on the statistics. 

The physics target, which was filled with liquid hydro- 
gen, was a 40-cm long cylinder with a radius of 2 cm. 
Continuous monitoring of the temperature and pressure 
permitted determination of the density with uncertainty 
of 0.2%. The target cell was surrounded by 24 "start 
counter" scintillators that were used in the event trigger. 

The CLAS detector utilized a non-uniform toroidal 
magnetic field of peak strength near 1.8 T in conjunction 
with drift chamber tracking to determine particle mo- 
menta. The detector was divided into six sectors, such 
that when viewed along the beam line it was six-fold sym- 
metric. Charged particles with laboratory polar angles in 
the range 8° — 140° could be tracked over approximately 
83% of the azimuthal angle. A set of 288 scintillators 
placed outside of the magnetic field region was used in 
the event trigger and during off-line analysis in order to 
determine time of flight (TOF) of charged particles. The 
momentum resolution of the detector was, on average, 
about 0.5%. Other components of the CLAS, such as 
the Cerenkov counters and the electromagnetic calorime- 
ters, were not used in this analysis. A more detailed 
description of the CLAS can be found in reference Q . 

The events were collected using a "two-track" trigger. 
This trigger required a coincidence between signals from 
the photon tagger and the CLAS. The signal from the 
tagger consisted of an OR of 40 of the 61 timing scintil- 
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lators, corresponding to photon energies above 1.58 GeV. 
This run setup where only part of the tagger was included 
in the trigger was intended to mainly detect photons 
above 1.58 GeV, avoiding the large bremsstrahlung con- 
tribution from photons of lower energy. However, events 
originating from a photon of energy below 1.58 GeV (cor- 
responding to a hit in counters 41-61) were recorded 
when they had an accidental (random) hit in one of the 
triggered tagger counters. By using these "accidental" 
events, we were able to collect data for photons with en- 
ergies from ~ 1 GeV to 1.58 GeV. This required applying 
an appropriate renormalization based on the probability 
for such a random coincidence to occur (see Sec. IVI|) . 
The signal from the CLAS required at least two sector- 
based signals. These signals consisted of an OR of any 
of the 4 start counter scintillators in coincidence with an 
OR of any of the 48 time-of-flight scintillators in the sec- 
tor. The rate at which hadronic events were accumulated 
was about 5 kHz; however, only a small fraction of these 
events contained the reactions of interest to the analysis 
presented here. 



III. DATA AND EVENT SELECTION 

The data reported here were obtained in the summer 
of 2004 during the CLAS "glla" data taking period, in 
which approximately 20 billion triggers were recorded. 
The relatively loose electronic trigger led to accumula- 
tion of data for a number of photoproduction reactions. 
During offline calibration, the timing of the photon tag- 
ger, the start counter and the time-of-flight elements were 
aligned with each other. Calibrations were also made for 
the drift times of each of the drift chamber packages and 
the pulse heights of each of the time-of-flight counters. 
Finally, processing of the raw data was performed in or- 
der to reconstruct tracks in the drift chambers and match 
them with time-of-flight counter hits. 

The reconstructed tracks were corrected for small im- 
perfections in the magnetic field map and drift cham- 
ber alignment, along with their mean energy losses as 
they passed through the target, the beam pipe, the start 
counter and air. In addition, small corrections were made 
to the incident photon energies to account for mechanical 
sag in the tagger hodoscope. 

The CLAS was optimized for detection of charged par- 
ticles; thus, the 7r + 7r _ 7r° decay of the r\ and the 7r + 7r~77 
decay of the rj were used to select the reactions of inter- 
est in this analysis. Detection of two positively charged 
particles and one negatively charged particle was re- 
quired. A 1-constraint kinematic fit to the hypothesis 
7P — > pir + ir~ (tt°/t]) was performed. This fit adjusts the 
momenta of all measured particles within their measure- 
ment errors such that energy and momentum are con- 
served and the missing mass is that of cither a tt° or 
an r\. The shifts in the momentum, combined with the 
known errors, yields a x 2 which is then converted to a 
probability (confidence level) that the event is the desired 



reaction. A cut was placed on the resulting confidence 
levels in order to select events consistent with one of the 
two topologies. Fits were run for each of the possible 
p, 7r + particle identification assignments using each of the 
recorded photons in the event. Photon-particle combina- 
tions with confidence levels greater than 10% were re- 
tained for further analysis. The trial identity as a proton 
or a pion for positive particles (assigned by the kinematic 
fit) was then checked using time-of-flight and momentum 
measurements . 

The covariance matrix of the measured momenta was 
studied using four-constraint kinematic fits performed on 
the exclusive reaction 7p — » pn + Tr~ in both real and 
Monte Carlo data samples. The confidence levels in all 
kinematic regions were found to be sufficiently flat and 
the pull-distributions (stretch functions) were Gaussians 
centered at zero with a = 1 (see Fig.[T|). The uncertainty 
in the extracted yields due to differences in signal lost 
because of this confidence-level cut in real as compared 
to Monte Carlo data is estimated to be about 3%. 

All negatively charged tracks were assigned a tt~ iden- 
tification. For positively charged tracks, the trial identi- 
fication from the kinematic fit was checked using time-of- 
flight and momentum measurements. The tagger signal 
time, which was synchronized with the accelerator radio- 
frequency (RF) timing, is used to obtain the start time 
for the event by accounting for the photon time of flight 
from the tagger to the reaction vertex. The stop time 
for each track was obtained from the TOF scintillator 
element hit by the track. The difference between these 
two times was the measured time of flight, t meas . Track 
reconstruction through the CLAS magnetic field yielded 
both the momentum, p, of each track, along with the path 
length, L, from the reaction vertex to the time-of-flight 
counter hit by the track. The expected time of flight for 
a mass hypothesis, m, is then given by 

The difference in these two time-of-flight calculations, 
Ato/ = t meas — t exp , was used in order to separate pro- 
tons from pions and to remove events associated with 
incorrect photons. 

Figure [2] shows Ato/ for tracks passing the kinematic 
fit under the 7r + hypothesis vs. Ato/ for the track pass- 
ing the fit under the proton hypothesis. The region near 
(0, 0) contains events where both tracks are good matches 
to their respective particle identification hypotheses. The 
2 ns radio-frequency time structure of the accelerator is 
evident in the out-of-time event clusters. Events outside 
of the black lines, where neither hypothesis was met, were 
cut from our analysis. This cut was designed to remove 
a minimal amount of good events. The Feldman-Cousins 
method [9( was used to place an upper limit on the sig- 
nal lost at 1.3%. Any remaining accidental events fell 
into the broad background under the 77/77', and were re- 
jected during the signal-background separation stage of 
the analysis discussed in Section ITVl 
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FIG. 1: (a) The confidence levels resulting from four- 
constraint kinematic fits performed on a sample of events to 
the calibration hypothesis jp — > ptt + tt~ integrated over all 
kinematics. The "peak" near zero consists of events that do 
not match the hypothesis, along with poorly measured (due 
to multiple scattering, etc.) signal events. Agreement with 
the ideal (flat) distribution for signal events is very good, (b) 
Example pull-distribution for the momentum of the tt~ from 
the same kinematic fits as in (a). Only events with a confi- 
dence level larger than 1% are shown. The line represents a 
Gaussian fit to this distribution. For this event sample, the 
parameters obtained are fi = -0.029±0.001, a = 1.086±0.001 
(the uncertainties are purely statistical), which are in very 
good agreement with the ideal values fj, = 0, a = 1. Both (a) 
and (b) are good indicators that the CLAS error matrix is 
well understood. We note that the 10% confidence-level cut 
used in the analysis corresponds to the relatively flat region 
of the confidence- level plot. (This figure is reproduced from 
reference [12] )■ 
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FIG. 2: (Color Online) Ato/ 7r + (ns) vs. Ato/ P (ns): Parti- 
cle identification cut for a sample of events that pass a 10% 
confidence level cut when kinematically fit to the hypothesis 
7P — > pir + Tv~ (tt°). The black lines indicate the timing cuts. 
Note the logarithmic scale on the intensity axis. 

the event. The lower-energy photon causes the reaction 
IV ~^ p7r + n~. When the event is reconstructed using 
the higher energy photon, there appears to be excess en- 
ergy and missing momentum along the photon direction. 
In some cases, this can be mis-reconstructed as a miss- 
ing 7T° moving along the beam direction. Our cut at 
cos#£ m . > 0.99 removes these events. A more detailed 
description of the entire analysis procedure presented in 
this paper can be found in Ref. [ig |. 

The resulting data have been sorted into bins in W. 
For the 77, there are 76 bins from W of 1.570 GeV to 
2.840 GeV. From 1.570 GeV to 2.100 GeV, the bins are 
10 MeV wide. From W = 2.100 GeV to 2.360 GeV, the 
bins are 20 MeV wide. In the range from 2.360 GeV to 
2.640 GeV, the bins are 40 MeV wide. Finally, there is 
a 50 MeV wide bin from W = 2.680 GeV to 2.730 GeV 
and a bin from 2.750 GeV to 2.840 GeV. The 77' data 
are divided into 44 bins from W = 1.900 GeV up to 
2.840 GeV. The bin widths for the 77' are the same as 
those used for the 77. 
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Fiducial cuts were applied on the momenta and an- 
gles of the tracks in order to select events from the well- 
understood regions of the detector. Included in these 
cuts was the removal of 13 of the 288 time-of-flight el- 
ements due to poor performance. In addition, events 
where the missing ir° was moving along the beam line, 
cos 8^ m > 0.99, were cut in order to remove leakage from 
the "fp — > p7r + 7r~ reaction (because of the very forward 
angle, the center-of-mass and lab angles are very similar, 
so a 0.99 cut in the center of mass corresponds to an even 
tighter cut in the lab frame). These arise from events in 
which two in-time photons occur where the higher-energy 
one is incorrectly associated with the charged tracks in 



IV. SIGNAL-BACKGROUND SEPARATION 

In order to extract physical observables from the 77 and 
77' photoproduction reactions studied here, it is neces- 
sary to separate background events from signal events. 
While such separation could be done via a simple side- 
band subtraction, our work on other reaction channels 
required the development of a more sophisticated event- 
based procedure. This event-based method, described in 
detail elsewhere [To|, [H| , was used to separate signal and 
background events. 

The procedure takes advantage of the fact that in 
the invariant mass distributions, the signal is a narrow 
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FIG. 3: Example bins demonstrating signal-background sep- 
aration. The plots are from the W = 2.11 GeV, cos#c.m. = 
0.75 kinematic bin for (a) r\ and (b) 77' . The unshaded his- 
tograms are all of the data in each 4-MeV wide bin, while 
the shaded histograms are these same events weighted by the 
background factors, 1 — Q. See text for details. 



structure, while the background is relatively featureless. 
The key feature of the procedure involves selecting each 
event's n c "nearest neighbor" events (we chose n c = 500) 

using the quantity A cos #c.m. as a distance metric be- 
tween events. Each subset of n c events occupies a very 
small region of phase space. Thus, the M 3w (M W7r ^) dis- 
tributions could safely be used to determine the proba- 
bility of each event being a signal event — the event's Q- 
factor. The M^j 'M wwrt distributions in each event's near- 
est neighbor event samples are fit to a narrow Gaussian 
(signal) plus a broad Gaussian and linear (background) 
function to determine the Q-factors. A few example in- 
variant mass distributions are shown in Fig. [3l Note that 
such a distribution and fit are generated for each event 
in order to determine the Q-factor for the event. These 



Q-factors are then used to weight each event's contribu- 
tion to the fits that are in turn used to determine the de- 
tector acceptance. The Q-factors are also used to weight 
each event's contribution to the differential cross section. 

Systematic studies were performed using different 
parametrizations of the background (including up to 
fourth order polynomials). From these studies, the sys- 
tematic uncertainty in the yield extraction due to the 
choice of background shape is estimated to be 4.1%(3.1%) 
for the 77(77') analysis. The point-to-point uncertainties 
obtained from the individual fits varied depending on 
kinematics; however, it was typically about 4%-5%. 



V. ACCEPTANCE 

The typical method of computing a detector accep- 
tance in a multi-particle final state such as those stud- 
ied here is an iterative procedure to achieve a physics 
model. To determine if the model is good, one generates 
a Monte Carlo sample thrown according to the model. 
These events are then passed through an accurate detec- 
tor simulation and the resulting event distributions are 
then compared to the physics distributions for the signal 
events as measured in the detector. After many itera- 
tions of the physics model, the two distributions should 
agree. At that point, one is able to use the physics model 
and the Monte Carlo sample to compute the detector ac- 
ceptance. 

There is an added complication to this procedure if 
there is background that is not easily separated from 
the signal. In such a case, one either needs to find a 
very clever way to effectively separate signal from back- 
ground, or include the background events in the physics 
model. Generally, one is required to iterate the physics 
and background models until satisfactory agreement has 
been reached. 

In this analysis, we have taken a more systematic ap- 
proach to computing the acceptance. First, we have used 
the Q-factors discussed earlier to allow us to produce dis- 
tributions of only signal events so that any model we use 
can ignore the background contribution. Second, we have 
used the fact that one can always expand any distribu- 
tion in terms of partial waves. Thus our procedure in- 
volves determining a partial wave expansion (our physics 
model), which will make the weighted Monte Carlo data 
agree with the observed signal sample. This partial wave 
expansion can be used to weight phase space generated 
events such that the Monte Carlo and data distributions 
agree. 

The efficiency of the detector was modeled using the 
standard CLAS GEANT-based simulation package and 
the Monte Carlo technique. A total of 100 million 77 
and 80 million 77' events were generated pseudo-randomly, 
sampled from a phase space distribution. Each particle 
was propagated from the event vertex through the CLAS 
resulting in a simulated set of detector signals for each 
track. The simulated events were then processed using 
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the same reconstruction software as the real data. In 
order to account for the event trigger used in this exper- 
iment (see Section [II]) , a study was performed to obtain 
the probability of a track satisfying the sector-based co- 
incidences required by the trigger as a function of kine- 
matics and struck detector elements. The average effect 
of this correction in our analysis, which requires three 
detected particles, is about 5%-6%. 

An additional momentum smearing algorithm was ap- 
plied in order to better match the resolution of the Monte 
Carlo data to that of the real data. Its effects were 
studied using four-constraint kinematic fits performed on 
simulated jp — > p7r + 7r~ events. After applying the mo- 
mentum smearing algorithm, the same covariance ma- 
trix used for the real data also produced flat confidence 
level distributions in all kinematic regions for the Monte 
Carlo data as well. The simulated 77 and rj events were 
then processed with the same analysis software as the 
real data, including the 1-constraint kinematic fits. At 
this stage, all detector and software efficiencies were ac- 
counted for. 

In order to evaluate the CLAS acceptance for the reac- 
tions 7p — ► pq and 7p — > prj' , we chose to follow the same 
procedure that we used in obtaining the acceptance for 
uj photoproduction [l^ ]. In this procedure, we expand 
the scattering amplitude for the pseudoscalar photopro- 
duction, Ai, in a very large basis of s-channel waves as 
follows: 

A^m^ ,m T ,m / (COS $^.m. ' ^0 
11 

«E E"<,^--/( cos ^'.U),(2) 

where a denotes a vector of 34 fit parameters, m^, m 7 , rrif 
are the target proton, incident photon and final proton 
spin projections on the incident photon direction in the 
center-of-mass frame, and A are the s-channel partial 
wave amplitudes. The s-channel structure of the am- 
plitudes, along with the details concerning the fit pa- 
rameters are described in Ref. [10|. The amplitudes are 
evaluated using the qft++ package [HI]. 

Unbinned maximum likelihood fits were performed in 
each W bin to obtain the estimators a for the parame- 
ters a in Eq. ^ . Background events were removed using 
the Q-factors directly in the fit as discussed in Refs. 10] 
and [111 ] . The results of these fits were used to obtain 
a physics weight, It, for each event. The weighted ac- 
cepted Monte Carlo events fully reproduce the real data. 
An example comparison is shown in Figure [4] for one W 
bin of the rj' data sample. The agreement between the 
weighted Monte Carlo and the data is very good. We 
note here that the results of these fits are not interpreted 
as physics, i.e. they are not considered evidence of reso- 
nance contributions to r)(r]') photoproduction. They are 
simply used to provide a complete description of the data. 

For a kinematic bin, the acceptance can be obtained 
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FIG. 4: (Color Online) An example fit result in a typical 
W bin for the r/ photoproduction. The data are shown as 
black squares, the phase-space accepted Monte Carlo events 
are shown as the blue dashed line, while the weighted Monte 
Carlo events discussed in Section [V] are shown as the solid 
red curve. The weighted Monte Carlo provides an excellent 
description of the data. 

as 

N acc 

acc{W, cos 9, h , f ) = , (3) 

j 

where N acc (N gen ) are the number of accepted (gener- 
ated) Monte Carlo events in the bin and the J's are the 
event weights discussed above. An accurate physics gen- 
erator would use the factors of / during the event genera- 
tion stage, rather than weighting the events. The result- 
ing acceptance calculation would be the same, modulo 
statistical fluctuations. 



VI. NORMALIZATION 

The measured rate of electrons detected by the tag- 
ger was used to compute the number of photons inci- 
dent on the target by sampling tagger hits not in coinci- 
dence with the CLAS. These rates were integrated over 
the live-time of the experiment to obtain the total pho- 
ton flux associated with each tagger element. Losses of 
photons in the beam line due to effects such as beam 
collimation were determined during dedicated runs us- 
ing a total-absorption counter placed downstream of the 
CLAS d|. 

The main method to calculate the experimental live- 
time during the "glla" run was based on a clock. This 
live-time calculation was checked by using the counts of a 
Faraday cup located downstream of the detector. While 
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the Faraday cup is a standard device for electron beam in- 
tensity measurements, in this case it was counting many 
fewer events produced by interactions of the photon beam 
with the target. Consequently, the statistical error of this 
second live-time measurement was high. However, the 
Faraday-cup-based measurement allowed us to observe 
that at maximum electron beam current, the actual dead 
time was about a factor of two higher than that given by 
the clock-based measurement. The high statistical error 
of the Faraday-cup-based measurement led to an uncer- 
tainty in the absolute live-time measurement of about 
3%. 

As discussed in Section [ill tagger counters 1-40 were 
in the trigger, while counters 41-61 were not. In order to 
have an event originating from a photon in the "untrig- 
gered" part of the tagger, the detector needed to record 
a random hit in one of the "triggered" counters during 
the trigger time window. The probability of such an oc- 
currence can be calculated using the electron rates. For 
the data here, we found that the probability of this hap- 
pening is 46.7%. This factor is then used to scale the flux 
in the untriggered counters. 

Events in the W — 1.955 GeV bin span the boundary 
between tagger counters 40 and 41. Because events in 
this energy range arose from both triggered and untrig- 
gered counters, the flux in this bin was deemed unreli- 
able. Thus, we report no cross sections for this energy. 
In addition to the above bin, the electronics in one of the 
tagger channels was not working properly during the run. 
This led to inaccurate flux measurements in the energy 
range W = 2.73 — 2.75 GeV. Differential cross sections 
will not be reported at these energies as well. 



VII. SYSTEMATIC UNCERTAINTIES 

An overall acceptance uncertainty of 5%-7%, depend- 
ing on center-of-mass energy, was estimated for this anal- 
ysis in Ref . [ijj ■ This includes uncertainties due to parti- 
cle identification (1.3%) and kinematic fitting confidence 
level cuts (3%), along with a relative acceptance uncer- 
tainty estimated by studying the agreement of physical 
observables obtained from each of the independent CLAS 
sectors (4 — 6%). The uncertainty on the normalization 
calculation was also estimated in Ref. [12] and found to 
be 7.9%. This includes contributions from photon trans- 
mission efficiency and live-time calculations. 

The acceptance and normalization uncertainties dis- 
cussed above were then combined with contributions 
from target density and length (0.2%), along with 
branching fraction (0.4% for the T), 1.5% for the rf) to 
obtain a total uncertainty, excluding the point-to-point 
contributions from signal-background separation (from 
the fits), of about 9%-ll%. The additional 4.1% and 
3.1% global signal-background uncertainties discussed in 
Section IIVI must then be added in quadrature giving to- 
tal uncertainties of about 10%-12% for both the rj and 
r}' . These errors are summarized in Table HI 



Error 


f) 


n 


A TOF Cut (PID) 


1.3% 


1.3% 


Confidence Level Cut 


3% 


3% 


Relative Acceptance 


4%-6% 


4%-6% 


Normalization 


7.9% 


7.9% 


Target Length 


0.2% 


0.2% 


Branching Fraction 


0.4% 


1.5% 


Background Shape 


4.1% 


3.1% 


Total 


10%- 


-12% 



TABLE I: A summary of the systematic errors associated with 
these measurements. The total presents the range of values 
seen across all data bins. 

VIII. RESULTS 

For the differential cross section calculations, da/dfl, 
each center-of-mass energy bin was divided into 20 bins 
in cos Sc.m? of width 0.1; however, results could not be 
extracted in every bin due to limitations in the detec- 
tor acceptance. In total, 1082 r\ and 682 r( cross section 
points are reported. The centroid of each bin is reported 
as the mean of the range of the bin with non-zero ac- 
ceptance. The results are shown in Figs. [H [7] and 
[5] and are available online in electronic form [lj| . The 
error bars contain the uncertainties in the signal yield 
extraction (point-to-point signal-background separation 
uncertainties and statistical uncertainties in the number 
of signal events), along with statistical uncertainties from 
the Monte Carlo acceptance calculations. The global sys- 
tematic uncertainties, discussed in Section IVII1 are esti- 
mated to be between 10%-12%, depending on center-of- 
mass energy. 

Both data sets have a forward peak that becomes more 
prominent with increasing energy, most likely due to 
some i-channel exchange mechanism. Both data sets also 
have a backward peak at our highest energies, this could 
be indicative of u-channel contributions. There are other 
interesting features in the data sets which might be evi- 
dence for resonance production; however, a partial wave 
analysis would need to be performed to determine if this 
is the case. In the r\ photoproduction, the most promi- 
nent of these are a bump near cos 9 = 0.2 in the low- 
energy data and a dip and a shoulder near cos 8 = — 0.1 
in the 2.1 to 2.6 GeV region. For the rj there is a similar 
dip and shoulder. 

Figure [5] shows a comparison of our rj results to pre- 
viously published CLAS data [5j. The previous CLAS 
results were published using cos6^? m bins of width 0.2, 
as opposed to the 0.1 width bins used in our work. To 
make this comparison, we have merged our bins (using a 
weighted average) to obtain the same binning as the pre- 
vious CLAS results. The agreement is good with most 
data points agreeing within the error bars. 
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FIG. 5: ^(/Ltb/sr) vs. cos#J? m for the jp — > prj reaction. Note that the vertical axis is linear for W up to 1.965 GeV and 
logarithmic above that. 
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FIG. 6: gn(/ub/sr) vs. cos#J? m . for the 7p — > pij reaction. Note the logarithmic scale on the vertical axis. 
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Figures [7] andJS] also shows the results from Tabs. II 
and IV of Ref. [l|| (a relativistic meson-exchange model 
which includes various resonance contributions). Five 
versions of the model are presented in [r|; each fit to 
the previously published CLAS data The previous 
data were not able to distinguish between the five ver- 
sions of the model. Near threhold, the Tab. II results 
clearly provide a better description of our data than those 
of Tab. IV. Given the poor performance of the models 
at higher energies, i.e. energies greater than what the 
models were fit to (W > 2.25 GeV), physics claims can- 
not be made from these models at this time. Our data, 
however, provide additional angular and energy cover- 
age, along with increased precision, which should allow 
for more reliable extraction of resonance contributions 
to rj' photoproduction by providing more stringent con- 
straints on the models. 

Figure [TU] shows a comparison of our rj results to pub- 
lished results from CLAS Q and CB-ELSA 0. The 
agreement between the three sets of measurements is fair 
to good at lower energies; however, at higher energies 
and forward angles large discrepancies begin to develop. 
The highest energy results reported by CB-ELSA are ap- 
proximately four times larger than our measurements in 
the most forward angles. Even for cos#J.' m « 0, the CB- 
ELSA results are approximately two times larger than 
ours at these energies. The figure also shows a SAID 
fit [TtJ to the earlier data. Given the differences between 
our measurement and the earlier data, the poor agree- 
ment between the fit and our results is not surprising. 
Including our results into the fit will likely have a signif- 
icant impact on the extracted physics. 

This extremely large disagreement at higher energies 
and forward angles motivated us to extract the "fp — > pr/ 
cross sections using an alternate procedure. This was 
carried out for points in our W — 2.46 GeV bin (which 
overlaps the second highest CB-ELSA bin). In this al- 
ternate procedure, we required the detection of only the 
p7r + , while ignoring the 7r _ . The 77 signal was then iden- 
tified in the missing mass off the proton. This topology 
has both a very different acceptance and a significantly 
larger background than the one in which all three charged 
particles are detected; however, the results obtained for 
the two topologies were in excellent agreement (see Fig- 
ure QT]). This study, along with the fact that our results 
from this data set in several other channels (sec Ref. 12] 
and [l8j]) are in good agreement with the world's data, 
has led us to conclude that it is very unlikely that there is 
some unknown acceptance or normalization issue present 
in our data. 

We can offer no explanation as to why the CB-ELSA 
results differ so much from ours; however, the self- 
consistency of the results obtained from our data set, 



using two distinct topologies, along with the high level 
of agreement with the world's data of cross sections ex- 
tracted for other reactions from this same data set pro- 
vide a high level of confidence in the results presented in 
this paper. 

Ultimately, one expects that these data, combined with 
other measurements, will facilitate a large scale partial 
wave analysis that will be able to identify the baryon 
resonance contributions to these cross sections. While 
we did attempt to carry out single-channel partial wave 
analyses of both of these channels [l|| (similar to that in 
Ref. [13]), the limited number of observables prevented us 
from drawing clear conclusions. Together with new mea- 
surements involving polarized beams and targets, these 
results should enable a deeper understanding of the nu- 
cleon resonances in the future. 



IX. CONCLUSIONS 

In summary, experimental results for 77 and 77' photo- 
production from the proton have been presented in the 
energy regime from near threshold up to W = 2.84 GeV. 
A total of 1082 r\ and 682 7/ cross section points are 
reported. The 77' results are the most precise to-date 
and provide the largest energy and angular coverage. 
The 77 measurements extend the energy range of the 
world's large-angle results by approximately 300 MeV. 
Unfortunately, discrepancies exist between the 77 results 
presented here and those previously published by CB- 
ELSA at higher energies. We look forward to seeing 
the impact these new results will have on existing models 
of baryon photoproduction. 
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FIG. 10: (Color Online) ^(ub/sr) us. W (GeV) for the 7p -> prj reaction, in cos0?. m . bins, from CB-ELSA Q(red filled 
circles), previous CLAS data [IJ(open squares), GRAAL 3] (light-blue open triangles), LNS 0] (purple open diamonds) and this 
work (blue open circles). The centroid of each bin in cos#c. m . is labeled on each plot. The solid line is a SAID [13] fit to the 
earlier data. 
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FIG. 11: (Color Online) ^ (Mb/sr) vs. cos 0? m : for the reac- 
tion 7p — > p?7. Differential cross sections near W = 2.46 GeV 
from CB-ELSA [2] (red filled circles) and from CLAS. The 
CLAS glla results presented in this work, which required 
detection of ptv + tt~ (blue open circles), and the "two-track" 
results discussed in the text, which only required the pix + to 
be detected (black filled squares) are in excellent agreement. 
The background levels in the "two-track" topology were too 
high to permit a reliable extraction of the signal yield for 
—0.4 < cos #c.m. < —0.2. The "two-track" error bars are 
purely statistical and do not contain any systematic uncer- 
tainty estimates on the signal-background separation. 



